Integrated application of multiple indicators and geographic information system-based approaches for comprehensive assessment of environmental impacts of toxic metals-contaminated agricultural soils and vegetables.

Conventional monitoring and mapping approaches are laborious, expensive, and time-consuming because they need a large number of data and consequently extensive sampling and experimental operations. Therefore, due to the growing concern about the potential of contamination of soils and agricultural products with heavy metals (HMs), a field experiment was conducted on 77 farm lands in an area of 2300 ha in the southeast of Shiraz (Iran) to investigate the source of metal contamination in the soils and vegetables and to model spatial distribution of HMs (iron, Fe; manganese, Mn; copper, Cu; zinc, Zn; cadmium, Cd; nickel, Ni, and lead, Pb) over the region using geographic information system (GIS) and geostatistical (Ordinary Kriging, OK) approaches and compare the results with deterministic approaches (Inverse Distance Weighting, IDW with different weighting power). Furthermore, some ecological and health risks indices including Pollution index (PI), Nemerow integrated pollution index (NIPI), pollution load index (PLI), degree of contamination (Cdeg), modified contamination degree (mCd), PIaverage and PIvector for soil quality, multi-element contamination (MEC), the probability of toxicity (MERMQ), the potential ecological index (RI), total hazard index (THI) and total carcinogenic risk index (TCR) based on ingestion, inhalation, and dermal exposure pathways for adults and children respectively for analyzing the noncarcinogenic and carcinogenic risks were calculated. Experimental semivariogram of the mentioned HMs were calculated and theoretical models (i.e., exponential, spherical, Gaussian, and linear models) were fitted in order to model their spatial structures and to investigate the most representative models. Moreover, principal component analysis (PCA) and cluster analysis (CA) were used to identify sources of HMs in the soils. Results showed that IDW method was more efficient than the OK approach to estimate the properties and HMs contents in the soils and plants. The estimated daily intake of metals (DIM) values of Pb and Ni exceeded their safe limits. In addition, Cd was the main element responsible for ecological risk. The PIave and PIvector indices showed that soil quality in the study area is not suitable. According to mCd values, the soils classified as ultra-high contaminated for Cu and Cd, extremely high for Zn and Pb, very high, high, and very low degree of contamination for Ni, Mn, and Fe, respectively. 36, 60, and 4 % of the sampling sites had high, medium, and low risk levels with 49, 21, and 9 % probability of toxicity, respectively. The maximum health risk index (HRI) value of 20.42 with extremely high risk for children was obtained for Ni and the HI for adults and children were 0.22 and 1.55, respectively. The THI values of Pb and Cd were the highest compared to the other HMs studied, revealing a possible non-cancer risk in children associated with exposure to these metals. The routes of exposure with the greatest influence on the THI and TCR indices were in the order of ingestion > inhalation > dermal. Therefore, ingestion, as the main route of exposure, is the route of greatest contribution to health risks. PCA analysis revealed that Fe, Mn, Cu, and Ni may originate from natural sources, while Fe was appeared to be controlled by fertilizer, and Cu primarily coming from pesticide, while Cd and Pb were mainly associated with the anthropogenic contamination, atmospheric depositions, and terrific in the urban soils. While, Zn mainly originated from fertilization. Findings are vital for developing remediation approaches for controlling the contaminants distribution as well as for monitoring and mapping the quality and health of soil resources.

Digital mapping and spatial modeling of some soil physical and mechanical properties in a semi-arid region of Iran.

The soil's physical and mechanical (SPM) properties have significant impacts on soil processes, such as water flow, nutrient movement, aeration, microbial activity, erosion, and root growth. To digitally map some SPM properties at four global standard depths, three machine learning algorithms (MLA), namely, random forest, Cubist, and k-nearest neighbor, were employed. A total of 200-point observation was designed with the aim of a field survey across the Marvdasht Plain in Fars Province, Iran. After sampling from topsoil (0 to 30 cm) and subsoil depths (30 to 60 cm), the samples were transferred to the laboratory to determine the mean weight diameter (MWD) and geometric mean diameter (GMD) of aggregates in the laboratory. In addition, shear strength (SS) and penetration resistance (PR) were measured directly during the field survey. In parallel, 79 environmental factors were prepared from topographic and remote sensing data. Four soil variables were also included in the modeling process, as they were co-located with SPM properties based on expert opinion. For selecting the most influential covariates, the variance inflation factor (VIF) and Boruta methods were employed. Two covariate dataset scenarios were used to assess the impact of soil and environmental factors on the modeling of SPM properties including SPM and environmental covariates (scenario 1) and SPM, environmental covariates, and soil variables (scenario 2). From all covariates, nine soil and environmental factors were selected for modeling the SPM properties, of which four of them were the soil variables, three were related to remote sensing, and two factors had topographic sources. The results indicated that scenario 2 outperformed in all standard depths. The findings suggested that clay and SOM are key factors in predicting SPM, highlighting the importance of considering soil variables in addition to environmental covariates for enhancing the accuracy of machine learning prediction. The k-nearest neighbor algorithm was found to be highly effective in predicting SPM, while the random forest algorithm yielded the highest R2 value (0.92) for penetration resistance properties at 15-30 depth. Overall, the approach used in this research has the potential to be extended beyond the Marvdasht Plain of Fars Province, Iran, as well as to other regions worldwide with comparable soil-forming factors. Moreover, this study provides a valuable framework for the digital mapping of SPM properties, serving as a guide for future studies seeking to predict SPM properties. Globally, the output of this research has important significance for soil management and conservation efforts and can facilitate the development of sustainable agricultural practices.